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Abstract 

The d-dimensional generalization of the point canonical transforma- 
tion for a quantum particle endowed with a position-dependent mass in 
Schrodinger equation is described. Illustrative examples including; the 
harmonic oscillator, Coulomb, spiked harmonic, Kratzer, Morse oscilla- 
tor, Poschl- Teller and Hulthen potentials are used as reference potentials 
to obtain exact energy eigenvalues and eigenfunctions for target potentials 
at different position-dependent mass settings. 

PACS numbers: 03.65.Ge, 03.65.Fd,03.65.Ca 

1 Introduction 

A position-dependent effective mass associated with a quantum mechanical par- 
ticle in the Shrodinger equation have attracted intense research activities over 
the years [1- 14]. It constitutes an interesting and useful model for the study 
of many physical problems. In the energy density functional settings to the 
many-body problem [1], the non-local term of the associated potential can be 
often expressed as a position-dependence of an appropriate effective mass M (r). 
Such an effective mass concept is used, for example, in the determination of 
the electronic properties of the semiconductors [2] and quantum dots [3], in 
quantum liquids [4], in ^He clusters [5] and metal clusters [6]. Nevertheless, 
within the Bohmian approach to quantum theory, the possibility of deriving 
the Schrodinger equation of particles with position-dependent effective mass 
from the Riemannian metric structure is explored and discussed (cf., e.g. [7]). 
Full and partial revivals of a free wave-packet, with position-dependent effective 
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mass, inside an infinite potential well are studied and documented [8]. TV-fold 
supersymmetry with position-dependent mass was reported by Tanaka [8], etc. 

However, in the study of Hamiltonians for particles endowed with position- 
dependent mass, M (r) = m.om{r), problems of delicate nature erupt in the 
process. The momentum operator, for example, does not commute with m{r). 
The choice of the kinetic energy operator is not unique, hence given rise of a 
quantum mechanical problem of long standing known as ordering ambiguity. 
Comprehensive details on this issue can be found in the sample of references in 
[8]- 

The above have formed, by large, the manifestos/inspirations of the recent 
studies on the one-dimensional Schrodinger equation for a particle with position- 
dependent effective mass [7-12]. However, conceptual and fundamental under- 
standings of the quantum physical phenomena may only be enlightened by the 
exact solvability of the Schrodinger equation. Yet, such exact solutions form 
the road-map for improving numerical solutions to more complicated physical 
problems. 

On the other hand, it is concreted that exactly solvable problems fall within 
distinct classes of shape invariant potentials (cf., e.g., [12,13]). Each of which 
carries a representation of a dynamical group and can be mapped into one an- 
other by a point canonical transformation (PCT) (cf., e.g., Alhaidari in [11,12] 
and Junker in [14]). The Coulomb, the oscillator, and the s-state Morse prob- 
lems, for example, belong to the shape invariant potentials carrying a represen- 
tation of so (2, 1) Lie algebra. In short, a PCT maintains the canonical form of 
Schrodinger equation invariant. 

In the PCT settings, one needs the exact solution of a potential model in 
a class of shape invariant potentials to form the so-called reference potential. 
The reference potential along with its exact solution (i.e. eigenvalues and 
eigenfunctions) is then mapped into the so-called target potential , hence exact 
solution for the target potential is obtained. Such a recipe is not only bounded to 
exact solutions but it is also applicable to the quasi-exact and conditionally exact 
ones ( cf, e.g., a sample of references in [15] on the quasi-exact and conditionally 
exact solutions), a consensus that should remain beyond doubts as long as the 
canonical form remains invariant. 

For the sake of completeness, the PCT approach for a quantum particle with 
a position-dependent effective mass M (r) = ruo m (r), in Schrodinger equation, 
should be complemented by its rf-dimensional generalization. Where interdi- 
mcnsional degeneracies associated with the isomorphism between angular mo- 
mentum £ and dimensionality d are incorporated through the central repul- 
sive/attractive core £{£+!) /r^ — > £d (^^ + 1) /r"^ of the spherically symmetric 
effective potential V^ff (r) ^ £{£ + 1) /r^ + V (r) (cf, e.g., [16-18] for more de- 
tails). To the best of our knowledge, the only attempt was made by Gang [18] on 
an approximate series solutions of the rf-dimensional position-dependent mass 
in Schrodinger equation. 

The forthcoming sections are organized as follows. In section 2 we provide 
the d-dimensional generalization of the point canonical transformation. We 
discuss the consequences of a power-law radial mass in the same section. lUus- 
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trativc examples are given in section 3. These examples include; the harmonic 
oscillator, the Coulomb, a spiked harmonic, a Kratzer-molecular, and a Morse 
oscillator as reference potentials. We also give two illustrative examples on 
the generalized Poschl- Teller and Hulthcn as reference potentials with different 
position-dependent singular masses in section 3. Our concluding remarks are 
given in section 4. 

2 PCT c?-dimensional generalization 

Following the symmetry ordering recipe of the momentum and position-dependent 

effective mass {M{f) = TOq m (r) , and a = 7 = 0, and /3 = —1 in equation (1.1) 
of Tanaka in [8]), the Schrodinger Hamiltonian with a potential field V (f) would 
read (in atomic units h = nio = 1) 

and assuming the (Z-dimensional spherical symmetric recipe (cf, e.g., Nieto in 
[18] for further comprehensive details on this issue), with 

* (r) = r-^'^-'y^Rn^,,, (r) Ye,,m, (0, , (2) 

Hamiltonian (1) would result in the following time-independent d-dimensional 
radial Schrodinger equation 

(3) 

Where = i+{d — 3) /2 for d > 2, i is the regular angular momentum quantum 
number, = 0, 1, 2, ■ ■ ■ is the radial quantum number, and m' (r) = dm, (r) /dr. 
Moreover, the d = 1 can be obtained through £4 = —1 and £4 = for even and 
odd parity, V = (—1)^''^^, respectively [17]. Nevertheless, the inter-dimensional 
degeneracies associated with the isomorphism between angular momentum £ and 
dimensionality d builds up the ladder of excited states for any given rir and 
nonzero £ from the £ = result, with that rir, by the transcription d^ d + 2£. 
That is, if En^^e {d) is the eigenvalue in d-dimensions then 

En^,e (2) = ^„.,^-i (4) = • • • = En^^i {2£) = K.,o (2^ + 2) (4) 

for even d, and 

£;„^,^ (3) = (5) = • • • = {2£ + 1) = £;„^,o {2£ + 3) (5) 

for odd d. Yet, a unique isomorphism exists between the S'-wavc {£ = 0) energy 
spectrum in 3D and in ID (i.e., En^^ (1) = En^fl (3)). For more details on 
inter-dimensional degeneracies the reader may refer to, e.g., [16-18]. 
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A substitution of the form R{r) = m{r)" (l){Z(r)) in (3) would result 
in Z' {r) = m{r) ~ manifested by the requirement of a vanishing coeffi- 
cient of the first-order derivative of (j){Z{r)) ( hence a one-dimensional form 
of Schrodinger equation is achieved), and Z' {r)^ = m{r) to avoid position- 
dependent energies. This, in turn, mandates v = 1/4 and suggests the following 
point canonical transformation 

q = Z{r) = f V^)dy =^ K,e. (r)) = m (r)"^/^ (r) . (6) 



Which in effect implies 



where 

t/,(,) = ^^-^^^ + "^^("^(^'/\ (8) 
8m (r) 32m (r) 4r m, (r) 

On the other hand, an exactly solvable (including conditionally-exactly or 
quasi-exactly solvable) rf-dimensional time-independent radial Schrodinger wave 
equation (with a constant mass ruo and cast in ?i = ruo = 1 units) 

+ ^iil^±ll + 2 [V (,) - .]} ^„^,, (g) = (9) 

would form a reference for the exact solvability of the target equation (7) . That 
is, if the exact/conditionally-exact/quasi-exact solution (analytical/numerical) 
of (9) is known one can construct the exact/conditionally-exact /quasi-exact 
solution of (7) through the relation 

Mi±i)+^(..)-a.M-E - + (10, 

Where Cd is the d-dimensional angular momentum quantum number of the 
reference Schrodinger equation. 

2.1 Consequences of a power- law mass m (r) = ar"^ 

With the radial position-dependent mass m (r) = ar'', the PCT function in (6) 
implies 

Z{r) = V^fy''/'dy=^^^r(^+'y' =^ t+^Z{r)=rV^) (H) 
and (8) gives 

1 (i%±4^U4fifii±ii^^ (12, 



16 V 2r2m (r) J 4 \ 2 (7 -|- 2)^ Z"^ (r) 
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Relation (10) in effect reads, with q = Z (r), 
with 



2g2 



A = -^ + l7 + 2rV4-^d(^rf + l) + (7-l)' + 27(3-d) 
Obviously, Eqs.(ll), (13) and (14) suggest that 7 = — 2 is not allowed. 

2.2 Remedy at 7 = —2 in a power-law mass m (r) = ar'^ 
For the case where m (r) = ar"^ equation (6) implies 



(14) 



q = Z (r) = \fa. J t ^dt = Valnr, 



(15) 



and hence 



t^47 = -2) 



f/4r,7--2)- 



id {id + 1) 



2a 



{id 



d-l 



2a 



(16) 



Which would only add a constant to the left-hand-side of (10) to yield, with 
£(i = and/or = —1 (i.e., only s-states and/or d = 1 states are available 
from the right-hand-side of (10) ), 

V{r)-Ud{^ = -2)-E ^ V{q)-e. (17) 

3 Illustrative examples 

3.1 m (r) = ar'^ with 7 7^ —2 

3.1.1 The harmonic oscillator reference potential 

The harmonic oscillator, 

V{q) = \xW. 

as a reference potential, with the exact d-dimensional eigenenergies and wavefunction 
enr,c, = A' (2nr + £d + , (18) 

i^n.,c, (9) = iMf'^' exp L^^+y^ {X\') (19) 
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respectively, would imply a target potential 



yir) = ^ar^^'; ^ = J^y (20) 

with corresponding d-dimensional eigenenergies and wavefunctions 

(7 + 2)0; 

E^^,e,= ^^ \' (2n, + A + l), (21) 

i?„.,.(r)=^„.,,(Cr)(*-)^-^exp (-^) ((CO^^) , (22) 
where 



A = A + 1/2 = |7 + 21-1 ^Aid {id + 1) + (7 - 1)' + 27 (3 - d) (23) 

and 

C= [2aa;/(7 + 2)]'/(^+'^ (24) 

It should be noted that this results, at = 3, collapse into Alhaidari's ones in 
example 5(a) of his Appendix in [12], where our w equals Alhaidari's C. 

3.1.2 The Coulomb reference potential 
The Coulomb, 

V (q) = -A/q, 

as a reference potential, with the exact d-dimensional eigenenergies and wavefunction 



^n.,c, (9) = N^rX. exp (^-h^^ (A„^,^^ ,) 



(26) 



respectively, would imply a target potential 

y(r)=--^r-i-^/2; C = A(7 + 2), (27) 
2va 

with corresponding rf-dimensional eigenenergies and wavefunctions 

_ gV2 1 
''"-^^-"(7 + 2)^n.+A+l/2)^ ^''^ 



ii„.,..(r)=A„,,,, (Cr)' exp -^^ ((^0 ) 




(29) 
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where 

ACT ] 1/(^+1) 

C = C K, 4) = (n. + A + 1/2)-^ 

(7 + 2) 



(30) 



It should be noted that our results (27)-(30), at d = 3, collapse into Alhaidari's 
ones reported in example 5(b) of his Appendix in [12], as his second solution 
of his Eq. (3.3) using the reference potential 3D harmonic oscillator. It seems 
that, in Alhaidari's second solution proposal of his Eq.(3.3) there is an implicit 
latent additional change of variables of a Liouvillean nature (cf., e.g., [24-26]) 
that led to a Coulomb-harmonic oscillator correspondence ( the reader may wish 
to investigate this issue following, e.g.,Znojil and Levai [26]). The proof of which 
is beyond our current proposal. 

3.1.3 A spiked harmonic oscillator reference potential 

A spiked harmonic oscillator ( or a Gold'man and Krivchenkov model), 

y(g) = AV/2 + /3g-V2, 

as a reference potential (cf,e.g.,[20]), with the exact d— dimensional eigenenergies 
and wavefunction 



3\ ;; 1 f„ 1^^ 



en^,c,=y [2nr + Cd + -y = - 3 + y (^^^d + - ) +/3. (31) 

i^n.,c, (a) = a„.,z:. (Ag)^^+^ exp Lif^'^ {X\^) (32) 

would lead to a target potential 

v r ^ 7+2 ^ ^ -7-2 2A^ - /3(7 + 2)' 



with corresponding rf-dimensional eigenenergies and wavefunctions 



(34) 



7-1-2' 



Rn.,. (r) = iV„.,, (Cr)(^+^)^+(-+^)/%xp (^-^J < {{Crr') , (35) 

It should be noted that Eq.(31) reduces to Eq.(12) of Yu and Dong [10] when 
d = 1, Cd = Q, —1, and 7 = — 3 to read 



En^fi = {2nr + 3) = V2A {2nr + 3) 
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where our = \f2A^ A is defined by Yu and Dong [10] as A = and our 

/3/2 = 15/8. One can also show that u in Yu and Dong is equal to [A'^g^] , and 
hence the corresponding wave function, for the d=\ case, 

(a) = Nn. exp (-|) Vl (u) , 

with A = 2, is exactly the same as that in Eq.(ll) of Yu and Dong in [10]. 
Moreover, our V (r) in (33) is the same as Eq. (4b) Yu and Dongs, of course 
with the proper amendments. 



3.1.4 A Kratzer's-type reference potential 

A Kratzer's-type molecular potential (cf, e.g., Fliigge in [21]), 

V{q) = -A/q + (3q-y2, 

as a reference potential with the exact rf— dimensional eigenenergies and wave- 
function 

-2 

^nr,Cd = "q ^rirXd = 7 : \"' (^6) 

« (nr+Cd + 1) 

i^n.,c. il) = (-^^) ^n^^' (v./:. q) (37) 

respectively, would imply a set of target potentials 

where C = ^(7 + 2), (3 = /? (7 + 2)^/4, and corresponding d-dimensional 
eigenenergies and wavefunctions 

(7 + 2) [nr + + 1/2) 



Rn.,i, (r) = An^,i, (^r)(*+^)^+'^ exp - 



(7? r) 



where 



i^(n. + .+ l/2)- 
(7 + 2) 



Lli{{vrr^'^') (40) 
]i/(i+i) 

(41) 
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3.2 m (r) = ar^ with 7 = — 2 

3.2.1 The spiked harmonic oscillator reference potential 

The spiked harmonic oscillator ( or a Gold'man and Krivchenkov model), 

V{q) = XV/2 + l3q-y2, 

as a reference potential (cf,e.g.,[20]), with the exact d— dimensional s-states' 
eigenenergies and wavefunctions 



3\ , 1 //l^^ 



^nr,o = X''[2nr + kd + ^j, kd = -^ + ^\^^j +/3. (42) 

would lead to a target potential 

y (r) = ^ {\nrf + (Inr)"' ; a = X-\ = ^ (44) 

2a 2 a 

with a corresponding d-dimensional eigenenergies and wavefunctions 

Er.^,,, = ^(2nr + fl+ ^^''^^^^^'^^^ j ,n = 2-Vl + 4aC2 (45) 



Rn. (r)=B„,-^(lnr)''+^/%xp 



(Inr)' 



Li 



((Inr)') , (46) 



It should be reported here that equations (45) and (46) reduce to the results 
obtained by Alhaidari (see example 5 in the Appendix of [12]) for id = and 
d = 3, However, it is worthy to mention that the ^d-dependence of the energy 
eigenvalues of the target potential are manifested by the consideration of the 
constant term in (16) of our proposal. 

3.2.2 A Kratzer's-type molecular reference potential 
A Kratzer's-type molecular potential, 

V{q) = -A/q + pq-y2, 

as a reference potential with the exact d-dimensional s-states' eigenenergies and 
wavefunctions 

-2 



(a) = N„.,o q"^^' exp ( ) Lf/+i (x„^,oq) (48) 
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would lead to the target potential 



V{r) 



A 



y/alnr 2a{lnrf 
with corresponding rf-dimensional eigenenergies and wavefunctions 



En 



{ed + kT + d-1 



2a 



(49) 



(50) 



Rnr (r)=Ar„,^(lnr)'=<*+' exp 



An,.o\/alnr 



xLf/+' (A„,,oV^lnr). 

3.2.3 A Morse-oscillator reference potential 
A Morse-oscillator potential of the form 

V {q) = yle-2°« - Be-"" ; B = 2A, 

as a reference potential (cf, e.g., [19]) with £4 = and/or £d = 
exact s-states rf-dimensional eigenenergies wavefunctions 



(51) 



^rir — ^ 



1 - 



1 



2Aa 



a = l/y/a 



-l.with the 



(52) 



u = \/8ale""«, s = ■yj-2aen^. 
would lead to a target potential 

nr)=-A(l-'^ 

with corresponding rf-dimensional eigenenergies and wavefunctions 



(53) 



2a 



1 / 1 

1 1= \'nr + - 

\f2AEi \ 2 



Rn. (r) =N^^ f - j e-«/2i^ (-n^, 2s + 1, «) . 



(54) 



(55) 



(56) 



Where u = V8aA/r. It should reported here that when d = 1 and dj, = 0,-1 
Eq.(55) reduces to Eq.s (14) and (!5) in [10]. 



10 



3.3 Two example on m (r) ^ ar'^ 

3.3.1 A Periodic Generalized Poschl- Teller reference potential and 

m (r) = a/4r {1 + r)^ q (r) = arctan-^/r 

A Generalized Poschl- Teller potential [21] of the form 



r (r — 1) >i{x — 1) 
cos^ C,q sin^ C,q 



(57) 



as a reference potential with the exact d-dimensional s-states' eigenvalues and 
eigenfunctions 

en^fi = \C[>^ + T + 2nrf (58) 



V'n, (?) =Cn^,^,T (sinCg)"^ (cosCg)^ 2-Fi {-nr,H+T + nr,>c+l/2]s,m^ Cq) 
would lead to a target potential 



V (r) = Vi{l+ r") + T/2 (^1 + - 
with d-dimensional s-states' eigenenergies and eigenfunctions 

Enr,td + T + 2nrf - (rji +V2 + %) 



(59) 
(60) 



(61) 



Where 



and 



a 



4r (1 + ry 



1 + r 



1 + r 



2F1 ( -nr,x + T + nr,>e+ 1/2; 



1 + r 



C = 1/Va, Vi = (24d - 9) /8a, r]^ = {8d - 9) /8a, 
V3 = - (32d - 22) /8a 



Vi = 



V2 = 



x(x-l)-(2a7?i+44 i£d + l)) 
2a 

T(T-l)-(2a7?2+44(^rf + l)) 

2a 



(62) 



(63) 



(64) 
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3.3.2 A Generalized Hulthen reference potential and m(r) = l/a^{r + 

1)2 =^ q{r) = a^^\n(r + 1) 

A Generalized Hulthen potential [22] of the form 

^(9) = -T^-7^ (65) 

as a reference potential with the exact d-dimensional eigenvalues and eigenfunc- 
tions (for the s -states) 

^nr,0 = 'y^I. (66) 



C„.,oe-Q-"^ £ (-1)^-^ ('^'^ + /3n^+ - - ^1 _ e-'^r 

(67) 

where Qn^ = 3(;ra — ^r)i and = l + 2Qn^ would leads to a target potential 
V{r) = ~ where a = a{a{d - l)/2 + 1) (68) 
with d-dimensional s -states' eigenenergies and eigenfunctions 

2 

En^fi = ^Qlr + ^ wliere B = a'^i'id - 3)/8 (69) 



Rn. (r) = C„.(l+r)-«"-V2 g - A + - - _ ,)'^ 

(70) 



4 Concluding Remarks 

In the point canonical transformation (PCT) method [23] ( an old Liouvil- 
lean change of variables spirit [24,25]) a Schrodinger-type equation often medi- 
ates ( via the existence of invertible parametrization of the real coordinates, 
r — > ''('z); and its few derivatives r' (q) , r" (g) , • • • ) a transition between 
two different effective potentials. In such settings, explicit correspondence (cf, 
e.g., Znojil and Levai [26]) between two bound state problems (i.e., the refer- 
ence/old and the target/new) is obtained. Within these Liouvillcan change of 
variables' spiritual lines, Alhaidari [11,12] has developed PCT-maps into tar- 
get/new position-dependent effective mass problems, in = 1 and d = 3. 

In this paper, a d-dimcnsional generalization of the PCT method for a quan- 
tum particle endowed with a position-dependent mass in Schrodinger equation is 
described. Our illustrative examples include; the harmonic oscillator. Coulomb, 
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spiked harmonic oscillator, Kratzcr-typc molecular, Morse oscillator, Poschl- 
Teller and Hulthen potentials as reference/old potentials to obtain exact energy 
eigenvalues and eigenfunctions for target/new potentials with different position- 
dependent effective mass settings. 

Finally, the applicability of the current PCT d-dimensional generalization 
extends beyond the attendant Hermiticity settings to, feasibly, cover not only 
PT-symmctric non-Hermitian Hamiltonians but also a broader class of r]-pscudo 
Hermitian non-Hermitian Hamiltonians [25-28]. This is already done in [29]. 
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